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Abstract. An iterative procedure for the explicit construction of the nontrivial subspace of all symmetry- 
adapted configurations with non-zero weight in the ground-state of the oo-dimensional Hubbard model 
is developed on the basis of a symmetrized representation of the transition operators on a sequence of 
Bethe-Lattices of finite depth. The relation ship between these operators and the well known mapping of 
the oo-dimensional Hubbard model onto an effective impurity problem coupled to a (self-consistent) bath 
on non-interacting electrons is given. As an application we calculate the properties of various Hubbard 
stars and give estimates for the half-filled Hubbard model with up to 0.1% accuracy. 

PACS. 7f.f0.-w Theories and models of many electron systems - 75.fO.Jm Quantized spin models 



1 Introduction 

The Hamiltonian of the half-filled Hubbard model is given 

by 



Bethe lattice this relation is given by [Q 



E 

i=i 
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(n-l) 
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The spin-dependence, a =] , |, has been explicitly retained 
in above formulas, in order include the antiferromagneti- 
cally ordered state. 



Once the parameters ej" l ' and l> have been ex- 
tracted from (H|) one needs to calculate (u>) by solving 
(pi) . The relevant Hilbert-space is of the order of 



(n-l) 



where Z is the coordination number of the lattice, i.e. 
Z = 2D on a simply hypercubic lattice. The model is well 
defined and nontrivial in the limit D — > oo Jl| . In this limit 
all correlations are local and the Hubbard model can 
be solved iteratively by mapping it to an effective Ander- 
son impurity model 0, which reads for the rfi 1 iteration 
as 



4 • 2 m t 
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where n = 0, 1,2 . . . denotes the number of iterations. The 
/| is the creation operator of the central (impurity) site 

and Mi™ -1 '' the number of state of the bath obtained 
from the previous (n — 1) iterations. The 'onsite' ener- 
gies e[™ and hybridization matrix elements V^™ ^ may 
be obtained from the local (impurity) Green's function of 
the previous iteration, G^™ (w). For the oo-dimensional 



The number of poles in G a (n)(u>) will of the same order 

of magnitude as (||), exponentially larger than Mj™ -1 ''. 
In any numerical treatment of it is therefore necessary to 
throw away an exponentially large part of the Hilbert- 
space JU|. 

Above route to solve the infinite-D Hubbard model is 
aimed to determine the one-particle Green's function and 
it is not obvious from (|f]) how to construct the eigenstates 
in the original Hilbert-space, in particular of the ground- 
state wavefunction. Here we show, that it is possible to 
construct the ground-state wavefunction iteratively, tak- 
ing the half-filled Bethe lattice as an example. 



2 Transition Operators on the Hubbard-Star 

A natural approach to solve (Q) iteratively is to consider 
the series of generalized Hubbard stars @] , which are trun- 
cated Bethe-lattices. We denote with S(0) the single site, 
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S(1) 





Fig. 1. Illustration of the Hubbard star S(l) and the star of 
the stars S(2). 



with S(l) the Hubbard star, with S(2) the star of stars, 
etc., see Fig. |l|. This approach has previously been ap- 
plied successfully to the study of the Mott-Hubbard tran- 
sition Q and allows to determine the ground-state of the 
infinite-D Hubbard model iteratively. In following we re- 
strict ourselves to the case of the antiferromagnetic state 
at half filling, generalization to other fillings and frustrated 
models are straightforward. 

We write the Hamiltonian of S(n) as 



+f/(/ T t / T -l/2)(/ i t / i -1/2), 



(4) 



where the /J are the electron creation operators at the 
central site and 



1 z 
At = — Vc* 



(5) 



Here the c\ a are the electron creation operators of the 

respective central sites of Sj(n-l). The in (^) are 

the Hamiltonians of the Sj(n-l), i = 1, . . . Z. 

We will show next that the knowledge of the exact 
eigenstates of H> n ~^ allows for the construction of the 
exact eigenstates of H' n > with a finite number of opera- 
tions, despite the fact that Z — > oo stars of order (n — 1) 
couple to the central site. 

We now introduce a notation for the eigenstates of 



H 



(n-l) 



Let \k m ,i) designate the k-th state in the sector 
m = (m+,m_) containing m^-particles of spin-a on sub- 
system i. We measure m a relative to the ground-state of 
m n-1 \ For the enumeration of the accessible states in the 
Hilbert space it is useful to define many-body transition 
operators 



Ifc 



m+s/ x^m 



(6) 



for the subsystem i with s = (1,0) or s = (0,1) re- 
spectively. We note that these operators obey the anti- 
commutation relations 



{ l^m+s')(^ml* ) |9m+s 

$i,3 3k' m ,q m+B \k m+s ,)(q m \i 
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Since A\ creates only symmetrized combinations of states 
on the periphery of the star, it is useful to rewrite A\ as 

Mr = E a (^m+s: ^m) l^m+sK^ml (8) 

where we introduced the symmetrized transition operator 

l^m+s)(^ml = ~~rr; l^m+s) (Ka\i (^) 

The coupling matrix elements in equation (^) are given as 

a (^m+sJ ^m) = (^m+sJ *I c 1,ctI*' ^m): (10) 

independent of i. 

Using Eq. ([?]) , commutation rules for the symmetrized 
operators are easily derived: 



( l^m+s) (^ml ) — r? l^m+s) (^mM^m+s) (K 



mlj 



= o. (11) 

Similarly we find: 

IOvWs'||?m+sK#ml = l fc m)( fc m+s'|ik m +s)(dj 
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Note that operators of type \ k' m ) (q' m \, which do not change 
the number of particles on periphery, are bosonic in na- 
ture, while the |fc m + s )(fcml are fermionic in nature, see 



3 Construction of the Hilbert Space 

We define the vacuum |0) of S(n) by 

|0) = |0)c®|0 (0l o),i), (13) 



where |0) c designates the vacuum of the central site and 
1 0(o.o) 7 l ne (niany-body) ground-state on the i-th leg of 
the peripheral sites. We define 

W) = fl\o), I U) = f]f\ |0) (14) 
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Fig. 2. Distribution of the \a(k± s , 0(o,o) )| 2 f° r S(2) as a func- 
tion of t*. Solid bars: s = (0, ±1), dashed bars: s = (±1,0). 
For U = 2t* +\a\ 2 are shown, for (7 = 8 — \a\ 2 are shown. The 



formation of the lower and upper Hubbard band for U 
can be seen. Here t* = 1 has been used. 



8t* 



on the same subsystem contribute to the ground-state of 

The other state, |fc(o,-i))(0(o,o)l I Tl)i m Eq. (|||) cou- 
ples to |0) and to the normalized state 

k(o,i)> (0(0,0)1 1^(0,-1)) (0(o,o) 1 1 t>- (19) 
The normalization of the state 

k(o,o)>(V-i)ll fc (o,-i)>( (o,o)ll T> = (20) 

k(o,o)> (0(o,o) 1 1 I) 



does again vanish for Z — > oo. 

We therefore conclude that of all states on the pe- 
riphery with m-f = = mi only the ground-state |0(o,o)) 
couples to the ground-state of H (™) . The complete Hilbert 
space is therefore spanned by 



n 



-(±1,0) 



>(°(o, 



{ k (±l,0)} 



II 

{ fc (0,±l)> 



where \x) = |0), |cr) or | fj.) and where all \k(±i 
{k(p,±i)} are mutually distinct, due to Eq. (pT|). 



(o,±i)) (0(o,o) (21) 
o)} and 



In order to iteratively construct the overall Hilbert 
space with nonzero coupling to the vacuum state defined 
above, we concentrate on the accessible states. Suppose 
for simplicity that the central site is occupied by an elec- 
tron with spin a = ]. The Hamiltonian if'") then couples 
the state | t) directly to the normalized states 

|fe(i,o))(0 (o ,o)l|0) and |A (0 ,-i)><0(o,o)| I U>- (15) 
The matrix elements of above states with | f ) are 

-t*o(ife(i,0)jO(o,o))/V2 and - t*a{k( 0i 



-i), 0(o,o) )/V2 



respectively, i.e. they are of order one in the limit Z — > oo. 

To construct the next set of states in the hierarchy 
of accessible configuration, we must act with A a on the 
periphery and re populate the central site. Without loss 
of generality we consider a = J.. For the first state in Eq. 
(jig), |fc(i,o))(0(o,o) I |0), we have only one possibility, 



k(o,-i))(0( O ,o)l 1^(1,0)) (O(o.o) I I I)) 
which is a normalized state. The normalization of 

k(i,-i))(fc(i,o)l |fc(i,o))(0(o,o)l I I) = 

1 =l<Z(i,-i))(0(o,o)IU) 



(16) 



(17) 



is 1/Z and therefore vanishes in the limit Z — * oo (here 
we have used (|l2|)). States with overturned spins on the 
same subsystem therefore do not couple to the ground- 
state of H^ n '. Generalizing this result one can show that 
only states on the periphery with 



|m T | + |m;| < 1 



(18) 



4 Mapping to the Anderson Model 

Let us define with Afi™ ^ the number of non-zero ma- 
trix elements a(k s , fc(o,o))- The Mi™ -1 ^ are related to the 
Mi n ~^ occurring in (|l]) via 



and 



M 1 ~ M (lfi) +M (-i,o) 



r (n-l) 



(22) 



(23) 



The hybridization matrix elements Vi of (|l|) are given 

by 



V, 



(n-l) 



1,1 



-t* 

72 



a(Z(±i,o)) 0(o,o))) 



(24) 



and respectively for a = [. The onsite energies e|™ of 
(PJ) are given by 



ff X) = E(l {±m )-E(0 {m ), (25) 



where the E(l s ) are the diagonal energies of H\ n X \ 
E{l s ) = (l s> i\H^\i,l s ), 



(26) 



independent of i. Note, that all energies are measured with 
respect to the Fermi-energy U/2, which is absorbed in Eq. 
(0) in the interaction term. 

Let's us see how things work out for S(l) and S(2). The 
half-filled ground state of S(0) is single occupied, let's say 
with an up-electron. We have then 

M (0) - - M (0) M (0) - 1 - M (0) 
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Table 1. The ground-state expectation values for the number 
of doubly occupied sites, (d), the local moment, (m) and of the 
kinetic energy, (T) (in units of t*), see Eq. (ps|). The values are 
for (7 = 8 and various Hubbard stars. Z is the weight of the 
states retained, see Eq. (^). L is the effective cluster size used. 





<d> 


(m) 


(T) 


Z 


L 


S(l) 


.007687 


-.492080 


-.249824 


1.0 


2 


S(2) 


.007797 


.491952 


-.249810 


1.0 


5 


S(3) 


.007807 


-.491941 


-.249896 


0.99993 


7 


S(4) 


.007807 


-.491940 


-.249901 


0.99995 
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with 

ej 0) = -U/2, 4 0) = U/2, V} 0) = -| = Vf > . 

The Anderson model for Hubbard star S(l) corresponds 
therefore to a two-site cluster || and the number of states 
contributing to the ground-state (which has one f- and 
one ^-electron) is Q (^) = 4. The numbers of one-particle 
and one- hole excited states for S(l) are all 

iW (±l,Q) _ Z ~ (0,±1) 

and S(2), the star of the stars, corresponds to a 5-site 
cluster. The ground-state is realized for three f- and two 
|-electrons (the state with two f- and three .[.-electrons is 
higher in energy). One has therefore 

M[% = 50 = <ar M S) = 100 = M (o%- 

S(3) corresponds therefore to a 151-site cluster. 

Above considerations are valid for constructing the 
exact ground-state and one-particle Green's function in 
the antiferromagnetic state. Effective Anderson models for 
S(l) and S(2) can although be constructed for the para- 
magnetic state, though with an increased size B . 



5 Discussion 

In view of the fact that the Hilbert-space increases ex- 
ponentially with every iteration one needs to discuss the 
feasibility of truncation schemes. From the anticommuta- 
tion rule 

C i,a C i,cr C i,a C i,a = 

for every i — 1, . . . Z and a =f , J, one can easily derive the 
sum-rule 

J2 Hh, (0 ,o))f + J2 H k ~ s > °(o,o))| 2 = 1. (27) 

which hold for both s = (1,0) and s = (0,1). Trunca- 
tion schemes become feasible, when an exponentially large 
number of the matrix-elements (|24|) become small due to 
the sum-rule (f27|). An indication of whether this is the 



Table 2. The same as in table Q for U = 4. 





(d) 


(m) 


(T> 


Z 


L 


S(l) 


.029127 


-.467646 


-.495066 


1.0 


2 


S(2) 


.030323 


.465713 


-.493618 


1.0 


5 


S(3) 


.030728 


-.465060 


-.495726 


0.99994 


15 



case or not may be seen by studying the distribution of 
the |a(/c ±s ,0 (o ,o))| 2 for S(2), which is given in Fig. |2|. 

Due to the formation of a local moment the matrix- 
elements shown in Fig. || are different for a = j. (solid bars) 
and a =t (dashed bars). Note, that Fig. || can also be 
interpreted as the one-particle Green's function of S(2) [Q. 
We observe that only a limited number of the transition 
matrix elements shown in Fig. ^|has an appreciable weight. 
The number of relevant poles increases with cluster size. 

We now consider in detail the ground-state properties 
of various Hubbard stars for U — 8, 4 and U = 2 (see 
Table @, § and |). For S(l) and S(2) we have constructed 
the exact ground-state wavefunctions and calculated the 
ground-state expectation values of the doubly-occupancy 
(d), the local moment, (to) and of the kinetic energy, (T), 
with (see Eq. (@) 

rf = /T/ T /l% ^=(/ T + / T -/J/ i )/2 

T =^ E (28) 
VZ <t=t,i 

Since S(3) corresponds to a cluster with L = 151 sites we 
cannot diagonalize S(3) exactly. For U = 2 and U = 4 we 
have approximated S(3) with a L = 14+ 1 site cluster, i.e. 
we have retained only the largest 14 transition matrix ele- 
ments (per spin) a(fc s , 0(0,0)) m Eq. (^) (all other matrix 
elements are set to zero). We also give in Table |^ and |^ 
the accuracy of this approximation, i.e. the contribution 
Z of the largest 14 matrix elements per spin to the sum- 
rule Eq. (H?]) . We observe that the truncation is better for 
larger values of U. For U = 8 it is possible to approximate 
S(3) already with by a L = 6 + 1 site cluster and S(4) by 
a L = 14+1 site cluster, see Table |. 

By setting the transition matrix elements to certain 
excited states on the periphery to zero a variational ap- 
proximation is obtained. It turns out that all the neglected 
excitations (those with very small a(k s , 0(o,o))) are high 
in energy in the sense that they are either above the up- 
per Hubbard band or below the lower Hubbard band, 
they contribute only to the tails of the respective Hub- 
bard bands. One can therefore estimate the contribution 
of these excitations to the ground-state of S(3) by per- 
turbation theory, their weight is < 10~ 3 for U = 2 and 
< 10~ 4 for U = 4, 8. Inspecting the data presented in the 
Table [ij ^ and || for various S(n) one sees that convergence 
with n is good, especially for [7 = 4 and U = 8. For U = 8 
the data has converged to within 0.01% for (to), to within 
0.1% for (T) and to within 1% for (d). 

In summary we have shown how the finite dimensional, 
minimal interacting Hilbert-Space of the Hubbard model 
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Table 3. The same as in table [j] for (7 = 2. 





<d> 


<m> 


(T) 


Z 


L 


S(l) 


.091221 


-.379144 


-.905348 


1.0 


2 


S(2) 


.092159 


.372893 


-.855795 


1.0 


5 


S(3) 


.097549 


-.359257 


-.878246 


0.99953 
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on the infinite-D Bcthe-lattice can be constructed itera- 
tively on the basis of symmetry adapted transition oper- 
ators. We explicitly provided the close relationship to the 
recursive construction of the one-particle Green's function 
based on the mapping to a self-consistent Anderson model. 
We have discussed the feasibility of truncation schemes in 
the iteration process, which are necessary due to the ex- 
ponential increase of the Hilbert-space at every step of 
the iteration. We have found that the sum-rule for the 
transition-matrix elements leads to a natural truncation 
criterion. We have applied the truncation scheme to vari- 
ous Hubbard stars and estimated the ground-state proper- 
ties of the half-filled Hubbard model up to 0.1% accuracy. 
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